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Significant  progress  was  made  in  the  application  of  PDF-based 
turbulent  combustion  models  to  the  Sandia  pilot-jet  nonpremixed  flames. 

These  PDF  calculations  were  described  in  two  recent  papers  (Xu  &  Pope 

2000  and  Tang,  Xu  &  Pope  2000)  and  they  are  summarized  in  this  report.  The  first 

of  these  works  showed  that  the  PDF  method  is  capable  of  describing, 

quantitatively,  the  phenomena  of  local  extinction  and  reignition 

(as  functions  of  jet  velocity  and  axial  distance)  that  are  observed 

in  these  flames.  The  second  paper  extended  the  methodology  to  include  radiative  heat  loss, 
and  the  calculation  of  NO  x. 
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SUMMARY 

Significant  progress  was  made  in  the  application  of  PDF-based  turbulent  combustion  models 
to  the  Sandia  pilot-jet  nonpremixed  flames.  These  PDF  calculations  were  described  in  two 
recent  papers  (Xu  &  Pope  2000  and  Tang,  Xu  &  Pope  2000)  and  they  are  summarized  in 
this  report.  The  first  of  these  works  showed  that  the  PDF  method  is  capable  of  describing, 
quantitatively,  the  phenomena  of  local  extinction  and  reignition  (as  functions  of  jet  velocity 
and  axial  distance)  that  are  observed  in  these  flames.  The  second  paper  extended  the 
methodology  to  include  radiative  heat  loss,  and  the  calculation  of  NOx. 

PDF  MODEL 

The  model  is  based  on  the  transport  equation  for  the  joint  probability  density  function 
of  velocity,  turbulence  frequency  and  thermochemical  composition.  For  the  combustion 
of  methane  considered,  the  chemical  kinetics  are  represented  by  the  augmented  reduced 
mechanism  of  Sung,  Law  &  Ghen  (1998).  Without  NOx  chemistry  (ARM1)  there  are  16 
species;  with  NOx  there  are  19.  For  computational  tractability,  the  effects  of  reaction  are 
implemented  by  the  in  situ  adaptive  tabulation  (IS AT)  algorithm  (Pope  1997). 

LOCAL  EXTINCTION 

The  piloted  jet  nonpremixed  flames  measured  by  Barlow  &  Frank  (1998)  consist  of  a  central 
jet  of  methane,  a  substantial  annular  pilot,  and  coflowing  air.  In  the  three  flames  considered, 
D,  E  and  F,  the  jet  velocities  are  Uj  =  50,  74  and  99  m/s,  respectively.  The  amount  of 
local  extinction  increases  with  the  jet  velocity.  _ 

Figure  1  shows  scatter  plots  of  the  CO2  mass  fraction  against  mixture  fraction  for  Flame 
F  at  30  jet  radii  downstream.  In  the  measurements  (left-hand  figure),  each  point  represents 
a  laser  shot:  in  the  calculations  (right-hand  figure),  each  point  represents  a  particle  from  the 
particle/mesh  method  used  to  solve  the  PDF  equation.  The  upper  curve  is  the  composition 
of  a  mildly  strained  laminar  flame.  The  fact  that  the  bulk  of  the  points  lie  beneath  this  line 
is  an  indication  of  local  extinction.  The  lower  curve  is  the  conditional  mean  from  the  data. 
Good  agreement  between  the  calculations  and  the  measurements  may  be  observed.  ■ 

Xu  &  Pope  (2000)  define  a  “burning  index”  (BI)  as  a  measure  of  (the  lack  of)  local 
extinction.  To  explain  the  definition  of  the  burning  index  we  again  refer  to  Fig.  1. 


The  vertical  dashed  line  shows  a  mixture  fraction  range  around  strichiometric  (defined  by- 
Barlow  &  Prank  1998).  The  lower  solid  circle  is  the  mean  of  the  data  points  that  lie  in  this 
range;  whereas  the  upper  solid  circle  is  the  value  from  the  laminar  flame.  The  burning  index 
is  defined  as  the  ratio  of  these  two  values.  Consequently,  a  burning  index  of  1  corresponds 
(roughly)  to  complete  combustion,  and  0  corresponds  to  complete  extinction. 


Experiment  results 


PDF  Calculations 


Figure  1:  Scatter  plots  of  CO2  mass  fraction  against  mixture  fraction  for  flame  F  at  an  axial 
distance  of  15  jet  diameters.  Left-hand  plot,  experimental  data  of  Barlow  &  Prank  (1998); 
right-hand  plot,  PDF  calculations  of  Xu  &  Pope  (2000). 

Figure  2  compares  the  measured  and  calculated  burning  indexes  of  different  species  (and 
temperature)  for  all  three  flames  as  functions  of  axial  distance.  As  may  be  seen,  with  few 
exceptions,  the  PDF  calculations  accurately  describe  the  dependence  of  BI  on  both  the  jet 
velocity  and  on  the  axial  distance.  These  results  demonstrate  a  significant  advance  in  our 
abilities  to  calculate  substantial  finite-rate  chemistry-turbulence  interactions. 


MAJOR  AND  MINOR  SPECIES 

For  flame  D,  Fig.  3  shows  calculated  radial  profiles  of  the  mean  and  r.m.s.  of  the  NO  mass 
fraction  compared  to  the  experimental  data.  Good  agreement  may  be  observed  for  both 
quantities,  and  the  downstream  evolution  is  well  represented.  It  may  be  seen  that  radiative 
heat  loss  does  not  have  a  significant  impact  on  NO  (in  this  flame). 

Figure  4  shows  the  conditional  mean  (conditioned  on  mixture  fraction  f )  of  various  quan¬ 
tities.  These  conditional  means  provide  a  good  indication  of  the  accuracy  of  the  chemistry 
calculations.  As  may  be  seen,  the  agreement  between  the  calculations  and  the  experiments 
is  uniformly  satisfactory.  The  trends  of  iVO  with  jet  velocity  and  downstream  distance  are 
accurately  represented;  and  the  accuracy  of  the  CO  calculations  shows  that  ARM  remedies 
previously  observed  discrepancies. 


Figure  2:  Burning  indexes  for  flames  D ,  E  and  F.  Filled  symbol:  experiments;  lines  with 
empty  symbols,  PDF  calculations.  Circle  and  solid  line,  flame  F\  down-triangle  and  dashed 
line,  flame  E\  up-triangle  and  dashed-dotted  line,  flame  D.  (From  Xu  &  Pope  2000.) 


Figure  3:  Radial  profiles  of  the  mean  (upper  plots)  and  r.m.s.  (lower  plots)  of  NO  mass 
fraction  in  flame  D.  Symbols,  experimental  data  (Barlow  &  Frank  1998);  solid  and  dashed 
lines,  PDF  calculations  without  and  with  radiative  heat  loss.  (From  Tang,  Xu  &  Pope  2000.) 
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Figure  4:  Conditional  means  in  flames  D  and  F.  Lines  and  symbols  as  in  Fig.  3. 
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